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We introduce a new class of exactly solvable boson pairing models using the technique 
of Richardson and Gaudin. Analytical expressions for all energy eigenvalues and first 
■ few energy eigenstates are given. In addition, another solution to Gaudin's equation is 

' also mentioned. A relation with the Calogero-Sutherland model is suggested. 

1. Introduction 

The notion of pairing plays a central role in the BCS theory of superconductivity. 1 
In 1963, R.W. Richardson showed that exact energy eigenvalues and eigenstates 
of the BCS pairing Hamiltonian can be computed if one can solve a given set of 
highly coupled nonlinear equations. 2,3 The limit of these equations in which the BCS 
coupling constant is very large were also obtained by Gaudin in 1976. Gaudin's tool 
was the algebraic Bethe ansatz method. 4,5 Gaudin's Hamiltonians are the constants 
of motion of BCS Hamiltonian in the large coupling constant limit. 

In this paper, we present a new class of exactly solvable boson pairing models 
using the technique pioneered by Richardson and Gaudin. The technique is outlined 
in the next section and in Section 3 we introduce the new class of exactly solvable 
models together with all energy eigenvalues and first few energy eigenstates. 

Richardson-Gaudin technique is based on an equation of which three solutions 
are known. To each solution there is a corresponding set of mutually commuting 
Hamiltonians and a related Lie algebra. In Section 4 we present another solution 
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to Gaudin's equation. We identify the corresponding set of mutually commuting 
Hamiltonians together with the related Lie algebra. 

2. Gaudin Algebra 

Gaudin starts with the following operators 4 : 

N 3 

3 = 1 a=l 

Here tog are arbitrary complex numbers and tf satisfy the SU(1, 1) commutators: 
[tf,tj] = -25 ij t° j , [tltf]=±5 l3 tf. (2) 
Then he shows that if we require these Hamiltonians to commute 

[/4°Uf]=0 (3) 
the coefficients wfj must obey the following equation: 

w ij w ]k + w ji w ik ~ w ik w jk = °- ( 4 ) 
Assuming i) the coefficients are antisymmetric under the exchange of i and j 

w% + w% = (5) 

ii) the coefficient wfj can be expressed as a function of the difference between two 
real parameters Ui and Uj, and iii) the operator X^i^i commutes with hj for all 
j, Gaudin found three solutions to Eq. (4). These solutions are given by 

wf, = for a = 1,2,3, (6) 

J Ui - Uj 

<=P^tm Ui - Uj )} ^f = sinhlp{ l_ Uj) y (8) 

The operators which are obtained by substituting these solutions in Eq. (1) are 
referred to as rational, trigonometric and hyperbolic Gaudin magnet Hamiltonians, 
respectively. 5 ' 6 In the rest of this section we consider the rational Hamiltonian 

The associated algebra, rational Gaudin algebra is an infinite dimensional complex 
Lie algebra generated by a one-parameter set of operators S + (X), S~(X) and S°(X). 
Commutators for rational Gaudin algebra are given as follows: 

[g+(A),g-(M)] = -2 g ° (A ]~ 5 ° (M) (10) 

A — /I 
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[g°(A),g ± CM)]=± g±(A |" g±(/i) (11) 

A — fJ, 

[S o (X),S a ^)] = [S ± (X),S ± ^)} = (12) 

Commutators of these generators for A = \i are found by taking the limit of the 
right hand sides of the above commutators as A — > [i. 
Let us define the operator 

H(X) = S°(X)S°(X) 1 -S+(X)S-(X) 1 -S-(X)S + (X). (13) 

Using the commutators (10-12), it is easy to show that H(X) forms a one parameter 
family of mutually commuting operators: 

[H(\),H(ri]=0. (14) 

One can diagonalize these operators simultaneously, starting from a lowest weight 
vector and using S' + (A) as step operators. Lowest weight vector |0 > by satisfies 

S"(A)|0>=0, and 5°(A)|0 >= W(X)\0 > (15) 

for every A € C. Here VF(A) is a complex function. One gets 

H(X)\0 >= E o (X)\0 >, E (X) = W(X) 2 -W'(X). (16) 

We write the Bethe ansatz state 

>=S+(^)S+^ 2 )...S+(U)\0> (17) 

for n arbitrary complex numbers £i, . . . , £ n S C. This is an eigenvector of H(X) 
if these complex numbers satisfy the following set of n Bethe ansatz equations 5 ' 6 : 

n 

w (&*)= E r~ T for « = l»2,...,n. (18) 

If ■ ■ ■ , £n S C is a solution of above equations then (17) is an eigenvector of 

H(X) with the eigenvalue 

IW) -*«•> . (19) 
a=i A ~ ?Q 

There exists a realization of the rational Gaudin algebra in terms of the SU(1,1) 
generators of Eq. (2) , given by 

Here u\,U2, . . . ,uat are arbitrary real numbers which are all different from each 
other and N > 0. In this realization H(X) of Eq. (13) is given by 

H(X) = V ^f 3 - . (21) 

^ K-A)(n, -A) 
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if (A) has simple poles on the real axis. Residues of —H{\)/2 at the points A = Ui 
are the rational Gaudin magnet Hamiltonians given in Eq. (9). Eq. (14) implies 



[H(X),h ( f ) }=0, and [hf\hf]=0. 



(22) 



3. A Model of Interacting Bosons 

The rational Gaudin algebra (10-12) can be realized in terms of the boson operators: 



ni 

1 J-L 
S-(A) = -X> 

p=l 



xi — A 

bpbp 
x\ — A 



ri2 + t 



9 5Z" C 



#2 - A 
2^ a ^X 

a=l 



«2 



/3=1 



x 2 - A 



(23) 
(24) 
(25) 



Here iri and X2 are two real parameters which differ from each other and = 
~E? a — 1. There are two types of bosons named bp and a a for — 1,2, ... , m and 
a = 1, 2, . . . , n 2 . The operators bp, bp and a a , a* a create or annihilate these bosons: 

[bp, b\] = 5pj [a a ,al] = S aa . (26) 

If we choose n\ — 5, n 2 = 1, = 1, and, E a = —1 then above operators can 

be considered as a generalization of the sd- version of the interacting boson model. 

7,8,940 

The lowest weight state is the boson vacuum and corresponding W(X) is 

1 / m n 2 



W(X) = 



4 \xi — A x 2 

In this realization, H(X) defined in Eq. (13) is a boson pairing Hamiltonian. 

For simplicity we subtract the ground state energy and write down H(X) 
H(X)-E (X): 

H(X) 



+ 



X 



(27) 



N d + N 



x\ — A 

+W(X) 



+ 



-S+(X)S-(X). 



(28) 



xi — A X2 — A 

Here Nd and N s are total number operators for the d and s type bosons, i.e. 

N d = Y. d W and N s = sU. (29) 



The last term in (28) is a pairing operator which corresponds to the Pq pairing 
operator of the sd-boson model. It involves the terms d)^d)^d v d v and s^s^ ss as well 
as the cross terms d^d^ss and s^s^d^d^. Therefore H(X) allows s-pairs to annihilate 
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each other to form s-pairs or <i-pairs. Similarly c?-pairs can also annihilate each other 
to produce s or d pairs. 

Since H(X) = H(X) — Eq(X), energy eigenvalues of H(X) can obtained from 
those of H(X) by subtracting the ground state energy. Eigenvalues of H(\) were 
presented in the previous section. Gaudin's equation (19) tells us that the energy 
eigenvalues of H (A) are 

^.-^ azsi , (30) 

Here, £i, £2, ■ • ■ , £n S C are obtained by solving the Bethe ansatz Eqs. (18). 

Here we present a method to calculate these energy eigenvalues in an analytical 
way without solving the Bethe ansatz equations. Let us start with substituting 
W(X) given in (27) into the energy formula (30) to find 



a=l 



-Si S 2 



(x\ — X)(xi — £ a ) {x 2 ~ X)(X 2 - 



(31) 



We see that it is not necessary to compute all the unknowns £1, £2, • • • , £n S C one 
by one to evaluate E n (X). It is sufficient to compute the following sums: 

Si = V 5 2 = V -. (32) 

We compute S 1 ! and ^2 using the symmetries of the equations of Bethe ansatz. It 
is easy to show that Eqs. (18) imply 

J2w(t a ) = and ^^ W ^ a ) = ^-J- (33) 

a— 1 a— 1 

no matter what the form of W(X) is. Substituting W(X) given in (27) in above 
equalities one finds two equations in Si and S2. Solving these equations for Si and 
S2 and then substituting them in equation (31) one finds 10 



1 



(xi - X)(x 2 — X) 



II 



E nW = 7 TT7 TT n(n - 1) + (m + n 2 ) - . (34) 



This is an analytical expression for all energy eigenvalues of H(X). It is computed 
without solving the Bethe ansatz equations. To find the corresponding energy eigen- 
states, however, one needs to solve the Bethe ansatz equations because each one of 
the quantities £1, & > • • ■ > £n G C is needed in order to write down the Bethe ansatz 
state in equation (17). 

Here, we present a method to turn the problem of solving the Bethe ansatz 
equations for the rfi 1 excited state into a problem of finding the roots of an rfi 1 
order polynomial. We substitute (27) in the equations of Bethe ansatz and then 
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apply the change of variables £ a = x 2 + ( a (xi — x 2 ). Bcthe ansatz equations then 
assume the following form: 

Er-V + T^-f^- fOT L > 2 "• (^) 

Ca - C/3 Ca 1 - Ca 

It was shown by Stieltjes 11 that when Ca obey (35), the polynomial 

n 

satisfies the hypergeometric differential equation 

C(l - C)K(C) - yK(C) + n(n + - l)j>„(C) = 0. (37) 

This means that Ca for a = 1, 2, . . . , n are the roots of the following polynomial: 
V* (-")*(" + n i/ 2 + "2/2 - l) fc ^ , . 

w^w k c =a (38) 

Here (a) = (39) 
(a)„=a(a + l)---(a + n-l) for n = 1,2, ■■■ (40) 

The polynomial in (38) terminates when k = n because (— n) n +i = 0. Consequently, 
finding the solution of the equations of Bethe ansatz for the n th energy level is now 
reduced to the problem of finding the roots on an n th order polynomial. In principle, 
we can find the roots of polynomials analytically up to the forth order. Therefore, 
first four excited energy eigenstates can be computed analytically. Solutions for the 
higher order eigenstates can be performed using numerical techniques. 

4. A New Solution to Gaudin's Equation 

In this section we briefly mention another solution to Gaudin's equation given in (4). 
Instead of the constraint (i) given in Eq. (5) this new solution obeys w°j + = 2q 
together with the other two constraints (ii) and (iii). The solution is given by 12 

wfj = q coth.[q(v,i — Uj)] + q for a = 1,2, 3. (41) 

Here q is a complex parameter. In the limit where q — > this new solution ap- 
proaches to the rational solution of Gaudin given by (6). For this reason we will 
denote the corresponding Gaudin magnet operators by h\ 9 ^ . In other words 

AT 

h l 9) = Ete coth W u * _ u j - )] + ^ • "*V ( 42 ) 

3 = 1 

These operators mutually commute: [h!f \ h^] = 0. This Hamiltonian and its Bethe 
ansatz solution was discussed in Rcf. 12. 
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5. A Connection with Calogero-Sutherland Model 

Calogero model was first introduced in 1969 as a many-body system in one dimen- 
sion with inverse square two-body interactions. 13 A variant of Calogero model intro- 
duced by B. Sutherland in 1971 includes two-body inverse sine square potentials. 14 
Both models are prime examples of exactly solvable systems and have many interest- 
ing features. Originally both models involved scalar particles. But generalizations of 
these models to particles with spin degree of freedom which are called spin-Calogero 
and spin-Sutherland models also received increasing attention. 15 > 16 > 17 > 18 > 19 Spin- 
Calogero model Hamiltonian is given by 

(X t Xj ) 



i=l i.j=l 



This model describes non-relativistic particles which are free to move in one- 
dimension. Here Xi denote the positions, pi denote the momenta and masses are 
scaled to unity, t i for i = 1, 2, . . . , N represent independent angular momenta of 
the particles. 

Spin-Sutherland model, on the other hand, is given with the Hamiltonian 

Hiq) = ^X> 2 +<^ E -Ju tj w (44) 

2 1~[ tl^i smh m x i - x i)\ 

i<3 

Here q is a complex parameter such that in the limit where q — > 0, the later model 
approaches to the former. It is possible to view Sutherland's model as a Calogero 
model on a circle, i.e. a system of particles which are constrained to move on the 
circle and interacting through Calogero potential along the cord distance between 
them. If we adopt this view then q = iir/d where d is the circumference of the circle. 

One can write the spin-Calogero model Hamiltonian using the rational Gaudin 
magnet Hamiltonians given in Eq. (9) as 

1 N 

ff (0) = ^£(^ + i fM)- (45) 

Similarly, the spin-Sutherland model Hamiltonian can be written in the same form 

i N 

tf (9) = ^£(^+*fb*,/ih) (46) 

i=l 

but this time we using the operators h\ 9 ^ which are given by Eq. (42). 

This connection between Calogero-Sutherland model and the Gaudin algebras 
is closely related to the Lax formulation of Calogero-Sutherland model. To see this 
let us first substitute wjj — — wfj = x(ui — uj) in Eq. 4. Here x(£) is a function 
which obeys x(£) + x(— £) = 2q. If we then take the derivative of Eq. 4 with respect 
to Ui and substitute £ = uj — Ui and r\ = Ui — Uk we find 

x(0x'(v) - AQx(ri) = x(£ + V )W(v) - *'(£)]• 
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It is well known that when x(£) satisfies this equation, a many body system with 
the two-body potential v = —x' is a Lax type system. In other words, its equations 
of motion can be written in the form of a Lax pair (see Ref. 20 for a review). If 
we use the rational solution of Eq. 4 given in Eq. 6 then the two-body potential 
v = x' is equal to the Calogero potential. On the other hand, if one uses of the new 
solution given by Eq. 41 then v = —x' turn out to be the Sutherland potential. 

6. Conclusion 

In this paper we suggested a few possible directions to extend the method of 
Richardson and Gaudin. Results of Section 3 are very encouraging because we 
see that one can obtain analytical results without explicitly solving the equations 
of Bcthc ansatz. The method of transforming the equations of Bcthe ansatz into a 
problem of finding the roots of polynomials is also introduced. Finding the roots of 
polynomials is much more established problem and this trick may come handy in 
numerical calculations as well. 

In this paper we relaxed one of the conditions imposed by Gaudin on the solu- 
tions of Eq. (4). This way we obtained a new and physically interesting solution. 
A new set of mutually commuting Hamiltonians are thus identified. Evidently, 
these Hamiltonians can be related to spin-Sutherland model in a way which is par- 
allel to the relation between the rational Gaudin magnet Hamiltonians and the 
spin-Calogero model. Application of algebraic Bethe ansatz technique to h\ q ^ yields 
a new Gaudin type algebra. This algebra admits a one parameter family of Hamil- 
tonians H(X) in a way which is analogous to other Gaudin algebras. 
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